Magnetic-field dependence of the localization length in Anderson insulators 
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Using the conventional scaling approach as weU as the renormahzation group analysis in d = 2 + e 
dimensions, we calculate the localization length £,{B) in the presence of a magnetic field B. For the 
quasi ID case the results are consistent with a universal increase of ^{B) by a numerical factor when 
the magnetic field is in the range £ <^ l^ ^ C(0), ^ is the mean free path, £^ is the magnetic length 
yjhcJeB. However, for d > 2 where the magnetic field does cause delocalization there is no universal 
relation between ^(-B) and ^(0). The effect of spin-orbit interaction is briefly considered as well. 
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Different universality classes, introduced by Dyson [Q , 
became the framework for considering the role of sym- 
metry breaking (due to the presence of a magnetic field, 
or spin-orbit scattering) in transport properties of disor- 
dered electron systems. This role has been well under- 
stood in the weak-localization regime, and in the vicin- 
ity of the Anderson metal-insulator transition |^-|4| . On 
the other hand, the classification in terms of the uni- 
versality classes appears to be quite useful also for quasi- 
one-dimcnsional systems in the strongly localized regime. 
One of the intriguing predictions for this regime ||^-g] is 
the existence of the universal relation, ^^ = /3 x C/3=i, be- 
tween localization lengths in different universality classes 
characterized by the parameter /3. There is, however, 
some controversy: while the transfer-matrix approach 
based on the maximal entropy ansatz ||7|,g| predicts that 
/3 = 1, 2, or 4 for orthogonal, unitary, and symplectic 
ensembles respectively, in accordance with the Dyson's 
definition for the random-matrix theory |H| , the exact cal- 
culation within the nonlinear a model [^ shows that were 
spin-rotation symmetry broken, /3 would be equal 4 even 
if time-reversal symmetry were broken as well, although 
such an ensemble would behave as unitary (/3 = 2) in rela- 
tion to the level statistics [Q. Nevertheless, the existence 
of a universal relation for f has been reliably established 
for quasi-one-dimensional systems, including quantum- 
chaotic ones iQ. However, the extension of the transfer- 
matrix approach to a higher dimensionality, as in Ref. 
[g[, is not obviously justifiable. 

In the present paper we analyze the infiuence of a 
magnetic field B and spin-orbit scattering (at i? = 0) 
on the localization length ^ applying the standard one- 
parameter scaling and renormahzation group (RG) ap- 
proach [H . The advantage of such a simple approach is 
possibility to apply the very same scheme to any dimen- 
sionality, d. We show that, while being in a qualitative 
agreement with the above-mentioned universality for ^ in 
quasi- ID systems, the scaling predicts the absence of any 
universal relation for d > 2. 

The RG equation for the dimensionless conductance 
g{L) is 



dlng/dlnL^ 13(g). 



(1) 



Here g{L) is defined for a d-dimensional system as 
the conductance in units of adie^/T^K) with ad = 
25(i_i/(27r)'^, Sd being the area of the d-dimensional unit 
sphere. The /3 function depends on the time-reversal and 
spin-rotational symmetry of the problem. Let us denote 
the (3 functions at i3 = 0, at large B and in the presence 
of "strong" spin-orbit scattering (at B = Q) respectively 
as /3o, Pu and /3s. The subscripts o, u and s stand for 
orthogonal, unitary and symplectic symmetries, respec- 
tively. 

The question is when is the symmetry effectively bro- 
ken. The answer depends on which length scale, L, we 
are in. We start at the microscopic scale L ^ £, with £ be- 
ing the mean free path, and go up to the range L ^ ^. At 
d > 2, that includes both localized and extended phases. 
In the extended phase ^ has the meaning of the corre- 
lation length which is indistinguishable in the critical 
regime, L ^ ^, from the localization length in the in- 
sulating phase. At d < 2, localization sets in with L 
increasing up to and beyond ^. 

In a semiclassical approximation, typical electron tra- 
jectories have phases (2Tr(j)/(j)g) (with 0^ = hc/e being 
the flux quantum) which are negligible for small fiuxes, 
when L <^ £^ = ^TicjeB and the problem is effectively 
that of i? = 0. On the other hand, once Bl? > (/i^, the 
phases are significant and the time-reversal symmetry is 
effectively broken. As long as localization has not set 
in, most of the relevant trajectories cover the whole area 
of the system, and increasing B beyond £^ ^ L should 
not further change the symmetry of the problem. Simi- 
lar considerations may well be valid for general Feynman 
paths, beyond the semiclassical approximation[J. Thus 
we are led to the following assumption: The correct (j{g) 



^ We do not consider here the range of larger magnetic fields 
^H ~ ^> where this treatment breaks down as the classical 
trajectories bend and orbit shrinking starts. So we leave out 



(for g ^ 1, L ^ ^j can be approximated by (3o for L '^ £^ 
and by j3u for L > 4f 

Similarly, we assume that for weak spin-orbit scatter- 
ing at i? = 0, f3{g) is approximated by /3o for L < 4o, 
and by /3s for L > £so (^so is the length of spin-orbit 
scattering) . 

The above assumption can be easily used in the RG 
procedure to estimate the change of ^ as a function of 
B. Let us start at d < 2 for weak disorder, k^i 3> 1, 
where localization still exists for large L. We start at 
a length scale L^, satisfying £ < Lg <^ ^, where g ^ 1 
and integrate the RG equation until getting to the regime 
where g ~ 1. This regime corresponds to the scale L ^ ^, 
and this is how ^ is determined. 

We approximate (3{g) by the weak-localization expres- 
sions lll,llll: 



Po{g) = d-2-l/g, 
Pu{g)^d-2-l/2g\, 
(3s{g) = d-2 + l/2g. 



(2a) 
(2b) 
(2c) 



The exact matching of P{g) at L ~ 4f or L ~ iso yields 
only unimportant corrections. We do not consider the 
case where both time-invariance and spin-rotations sym- 
metries are broken. 

We begin with a quasi ID wire, having N^ conduction 
channels (i.e. N^ transfer states below Ep). For the wire 
of a 2D cross section the classical Drude conductance is 
given by 



g,,{L) = ^N^j 



(3) 



We start integration of Eq. ([I]) from a small L^ <C ^ 
(but Lq :^ i). For B large enough, ^^ < ^o (Co is the 
B = localization length), we first integrate Eq. (ll) in 
the orthogonal symmetry range, from L^ to 4r, with the 
/3-function /3o, Eq. (pa[), to obtain (using gci{Lg ^ 1) 



g{L) + l^g,t{L), 



L<L 



(4) 



Note that Lg has cancelled out. For L's above 4, , starting 
with L = ijj we integrate Eq. (|^) with the /3-function /3„, 



Eq. (^, to find (using 5(4) > 1) 

g^{L) + 1/2 = gl,{L), L>£„ 



(5) 



Note that 4f has cancelled out as well, this will be seen 
to be the source of the universal increase of ^ with B in 
the whole regime considered in the quasi ID case. 

For a very weak magnetic field, 4? ^ Co, Eq. (0) is valid 
in the whole range of L, up to Co, and the behavior is "or- 
thogonal" throughout, still allowing small i3-dependent 



of account interesting effects that may occur in the range of 
even higher fields, ^^ ^ £-/kp, where one has to scale g^x and 
Qxy together |ll. 



corrections to C- At L = Co, g{^o) = A, A being a con- 
stant of order unity. Thus 



Co 



TT NJ 

SA + l 



(6) 



in accordance with the well-known Thoulcss result [[l4[ . 
In the purely unitary case, i.e. for a "strong" magnetic 
fields 4f ^ Lg, one has an increase of C: 



C -- 



NJ 



3 ^v42 + 1/2 



(7) 



A similar increase takes place in the presence of spin-orbit 
scattering. Repeating the same procedure as above, one 
finds that for the purely symplectic case 



6 = ; 



NJ 
3A~l/2 



(8) 



While the precise value of C«/Co and Cs/Co cannot be de- 
termined by this simple approach, this increase is quali- 
tatively due to (3{g) being larger in the unitary case, and 
even more so in the symplectic case so that one has to go 
to larger scales to get g ~ 1. Since i^ dropped from Eq. 
(^), the increase from Co to Cm is universal (i.e. it does not 
depend on B) in the whole range £ '^ £„ ^ C- Similarly, 
the increase of C in the symplectic case does not depend 
on the spin-orbit scattering strength for £ ^ 4o ^ C- 
Thus our scaling consideration, while not able to yield 
the constants, is consistent with the exact result [g that 
in a quasi ID metallic wire there is an increase of C in 
the appropriate, rather wide, range of B. 

The physical reason for the increase of C with B is now 
very clear. It is due to the well-known delocalization for 
the weakly localized regime L <^ ^o, caused by the elim- 
ination of coherent time-reversed paths by the magnetic 
field B. The scaling procedure enables us to approxi- 
mately carry this effect all the way up to L ~ C and see 
by how much C is changing. 

The same procedure which leads to semi-quantitatively 
correct results in quasi ID and illuminates the physics of 
delocalization by the field, can be repeated for a 2D film. 
We start here at Lg ~ € with a conductance gg and obtain 
for Co > 4 ^ ^ 



C = 4exp(.g2) 



(9) 



where g^ = g{£„) = 5o ~ ln(4/^)- This ranges from the 

orthogonal value given by gg ~ ln(Co/^), i.e. Co — £e^° for 

2 
4f ^ Co, to the unitary one, Cm — ie^° for 4 ~ ^- Since 
gg 3> 1, Cm is overwhelmingly larger than Co- The ratio 
Cm /Co is obviously not universal - it changes by many 
orders of magnitude in the range ^ ^ 4 ^Co- 
in the symplectic case, the (3s function, Eq. ([2q), is 
always positive in the weak localization regime. So for 
£,0 < ^so, there is an "antilocalization" effect: the con- 
ductance is decreasing with increasing the scale up to 
L ^ 4o- Then the /3 function in Eq. (nh changes sign 



from negative, Eq. (2a), to positive, Eq. (|2q), and g be- 
gins to increase. The antilocalization character of the /3 
function might change at 5 ~ 1 due to further terms of 
the 1/g expansion |15t] . Such a change, if any, is definitely 
beyond the scope of the present treatment. 

We now consider the situation for rf > 2, where the 
Anderson transition exists at a certain g — gc- Here we 
use the RG equation (|^) hnearized near gc (for g—gc <C 1) 
with 



V gc 



(10) 



where v is the exponent of the diverging ^ near the tran- 
sition. Indeed, the solution to Eq. (|l|) with the linearized 
/?-function (nO[) for g {— g{L)) at a general scale L is 
written as: 



9c 



9= {9c-9o){L/L^f''' 



(11) 



where g„ = giL^). 
i9o < 9c) as 

g = gc\l- {L/0 



We rewrite this in the localized phase 



i/i, 



^ = -^0 (1 - 9o/9cY 



(12) 



From now on we take the above notation, with the con- 
stants gc and i/, to describe the orthogonal case {B = 0) 
with ^ — ^o- The unitary ("large i?") case satisfies the 
same scaling theory but with the constants gc and v re- 
placed by gcH and v^ . In the e expansion for d = 2+e, the 
transition at small e occurs in a range describable by per- 
turbation theory, and we can use the weak-localization 
expressions (0) to obtain the following constants 



9c 



1/e, ly = 1/e; .g,„ = (l/2e) 



1/2 



^„ = l/2e. (13) 



We shall assume that gcH < 9c also at e=l and use our 
ansatz that the system behaves according to orthogonal 
or unitary symmetry for L < i^ ot L > i^ respectively. 
We shall also assume that the crossover is describable by 
the linearized P{g). The latter assumption, too, is justi- 
fied in the e-expansion. We do the considerations below 
for a general d > 2. They are certainly valid for e ^ 1 
and should be qualitatively correct at d = 3. (Unfortu- 
nately, such an approach is not feasible for the spin-orbit 
case where gc is not large even for small e, and the e- 
expansion could hardly provide with useful insight into 
the 3D case, besides the fact, though, that the spin-orbit 
interaction still does cause delocalization for £ '^ £so '^ £, 
in the weakly localized phase near the transition.) 

Again, for ^ ^ 4r ^ ^1 ^'^ use the orthogonal results, 
Eqs. (|ll|), (111), from L„ = ^ to 4: 

5(4) =9c- {9c - 5o) (4/^)* = 9c [1 - (4/^0)*] (14) 

At the scale L ^ i^^, we use the unitary results. When 
5(4) > 9cH, the decrease in g, Eq. (M|), is changed by an 
increase as the unitary /^-function (00^ is positive at this 
scale. Then, the critical conductance as function of B, 



gc{B), is given by g^ for which q(C) = gen- Substituting 
into this equation g(4)) El- (|l4|), and solving for g^ = 
gc{B) we find: 



5c(i3)=5c-(<?c-5c«)(44)" 



(15) 



To check this, we see that gc{B) = gen for i^ = £, and 
gc{B) -^ gc for £„ — > 00 (i? — > 0). The real divergence of 
^ near the transition occurs for a given B when gdB) — > 



C« = 4 1 



9cH 



At this scale, g is given by Eq. ( [12|) with substituting 
9cH, £,HJ and v^ for respectively gc,£,o, and v. Clearly, 
there is nothing universal in the ratio of S.h/^- A sim- 
ilar conclusion for d ^ 2,3, both for the magnetic and 
spin-orbit case, has been made in Ref. |1^ based on the 
heuristic picture of localization of Ref. 1 17J . 

The most important aspect of delocalization by a mag- 
netic field in the range considered (^ < 4f ^ 'fo) at d > 2 
is the shift of the transition towards delocalization. This 
is seen from the fact that for small e, gen *S 9c and it is 
expected that gen < 9c also at d=3. For a given field B 
there are three regimes: for g < g^B), there is a localiza- 
tion also with the magnetic field but S, is larger (arbitrar- 
ily larger near gc{B)); for gdB) < g < gc the magnetic 
field destroys localization completely; for g > gc, there 
is no localization but the correlation length ^ appears to 
be larger at i? = (especially just above gc). The shift 
of the transition with field has been discussed in Refs. 
p8| , p9| . Clearly, the effect of B for d > 2 is not universal 
(except, of course, that the asymptotic critical behav- 
ior is governed by the unitary universality class for any 
B). It is only for quasi ID systems that ^ increases by 
a universal factor in the whole range ^ ^ 4 ~ Co- It 
is interesting that the considerations based just on the 
symmetry breaking within the transfer matrix approach 
W seem to be valid only for quasi ID systems. They had, 
however, highlighted the very interesting and at the time 
almost unexpected possibility of substantial delocaliza- 
tion by the magnetic field. 

The scaling cross-over considerations exemplified here 
can be applied to various physical situations involving 
the phase-breaking length £(j,, alongside with 4f and ^so in 
appropriate ranges of the localized phase. This produces 
a variety of interesting magnetoconductance phenomena 
in the localized phase. 
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